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We present numerical and analytical results for the reflection and transmission properties of matter 
wave solitons impinging on localized scattering potentials in one spatial dimension. Our mean field 
analysis identifies regimes where the solitons behave more like waves or more like particles as a result 
of the interplay between the dispersive wave propagation and the attractive interactions between 
the atoms. For a bright soliton propagating together with a dark soliton void in a two-species Bose- 
Einstein condensate of atoms with repulsive interactions, we find different reflection and transmission 
properties of the dark and the bright components. 
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I. INTRODUCTION 

Solitons are exceptionally stable wave phenomena that 
appear in a variety of physical systems. They propagate 
without spreading due to a fine balance between the dis- 
persion and the non-linearity of the physical system. Af- 
ter their discovery as stable waves in shallow waters by 
Russell in 1834, solitons have been found in very different 
non-linear systems. As an example, both bright optical 
solitons which are peaks in the intensity and dark soli- 
tons which are intensity minima have found applications 
in optical transmission through fibres. In this article we 
consider soliton waves in Bose-Einstein condensates of 
neutral atoms. 

The condensation of a cold gas of atoms permits a 
description of the many-body problem by a single con- 
densate wave function order parameter, which is in turn 
described by the non-linear Gross-Pitaevskii wave equa- 
tion [1^. The non-linearity in Bose-Einstein condensates 
arises from the repulsive or attractive interactions of the 
particles in the condensate. In equivalence with non- 
linear optics, dark and bright solitons are, respectively, 
localized dips and peaks in the Bose-Einstein condensate 
density. 

The dynamics of solitons in inhomogeneous conden- 
sates with slowly varying potentials have been investi- 
gated and is adequately described by a particle-like be- 
havior [2HS]- In this article, we shall further elaborate 
the analysis of matter wave solitons in ID-condensates, 
investigating their scattering properties on localized po- 
tentials. 

In Sec. II, we introduce the Gross-Pitaevskii equation 
and a convenient mapping to dimensionless coordinates 
that will be used throughout the paper. In Sec. Ill, we 
introduce the bright, the dark, and the composite dark- 
bright solitons, and we briefly review their properties. In 
Sec. IV, we show that scattering of bright solitons on a 
low and wide potential barrier can be well understood as 
the motion of a classical particle governed by an effective 
Lagrangian, extracted from the quantum problem, while 
the scattering on high and narrow barriers can be un- 
derstood from the single particle quantum problem. In 



Sec. V we present results for the scattering of a dark- 
bright soliton on a potential felt only by the bright soli- 
ton atoms, and we show to what extent their scattering 
is transferred to the dark soliton. Sec. VI concludes the 
paper. 



II. MEAN FIELD THEORY OF INTERACTING 
ATOMS 

We treat the Bose-Einstein condensate of atoms at zero 
temperature in a Hartree mean field theory, which ne- 
glects particles outside the condensate and treats short 
range interactions in the gas by an interaction potential 
which is proportional to the mean atomic density. This 
yields the ID timedependent Gross-Pitaesvkii equation 
(GPE) 

where the order parameter, the macroscopic wave func- 
tion, is normalized to the total number of atoms 

y dV \ij[r)\^ = N (2) 

so that \tp{x,t)\'^ represents the particle density. The ID 
interaction strength 17 is a known function of the 3D low 
energy scattering properties and the transverse confine- 
ment of the atoms [T]. 

In an attractive BEG with g < 0, (meta-)stable bright 
soliton solutions with a finite number of particles exist, 
while in a repulsive BEG with g > 0, dark solitons may 
exist in the form of holes in the background density. 

For numerical purposes it is advantageous to rewrite 
the GPE in a dimensionless form, and for a homogeneous 
condensate with repulsive interactions, the linear atomic 

density, no and the healing length ^ = \/ J^l„ provide 
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natural scale factors. We thus introduce dimensionless 
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quantities for length, energy and time, labelled by a tilde 

(3) 



and obtain the equation 
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where the dimensionless wave function 
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1 in the ho- 



is normalized to unit bulk density 
mogeneous part of the condensate. 

A condensate with attractive interactions has no nat- 
ural length scale like the healing length, but for conve- 
nience in the following, we introduce an arbitrary length 
unit which may be chosen to fit a relevant experimental 
length scale, and we set 
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t = t- 
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Inserting this in the GPE, Eq. ([T]), and introducing the 

\g\ ^ |V'I^ leads to the 



rescaled wave function, 
dimensionless GPE, 
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The rescaled number of particles in the condensate is 
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Thus for given N, the choice of N determines the length 
scale I and sets the dimensions of the problem. 

Due to the presence of the non-linear term it is in 
general necessary to solve the Gross-Pitaesvkii equation 
by numerical means. The numerical results presented 
in the following are based on propagation of the one- 
dimensional, dimensionless and non-homogeneous GPE 
using the Crank-Nicolson method. For brevity we will 
omit the tilde above dimensionless quantities through- 
out. 



III. SOLITONS 

In this section we recall the soliton solutions of the 
homogeneous Gross-Piaevskii equation. 



A. Dark solitons 

Although the scattering properties of dark solitons are 
not discussed in this paper we will first introduce the 
dark soliton. This will serve as a background for the 
later discussion of the dark-bright soliton. 

Dark solitons are formed in repulsive BEGs, governed 
by the homogeneous GPE, Eq. (|4|, and fulfilhng the 
boundary condition lima;_^±oo |'0(2;)| = tiq, where no is 
the bulk value of the condensate density and also the 
dimensionless chemical potential. A dark soliton soliton 
moving with velocity v is then described by 
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where = ^/y^l — with ^ the healing length and 

c — the speed of sound, both unity in our dimen- 

sionless units [1]. The dark soliton has a density dip 
around x = vt, and we will characterize the size of the 
soliton by the number of atoms missing, Nd, compared 
to the homogeneous case. 



iVo = Jdx (no - |i^d|') = 2710^"^ 1 - 



(10) 



The density depression moves in the opposite direction 
of the atoms in the condensate, whose motion is given by 
the gradient of the phase of the condensate wave function. 
The energy of the dark soliton is given by the difference 
between the energy of the condensate containing a dark 
soliton and the energy of the homogeneous condensate. 
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The soliton energy decreases with increasing soliton ve- 
locity, and for ^ <C 1 one finds 
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corresponding to an effective negative mass of —Auq^. 



B. Bright solitons 

If the interactions in the BEG are attractive {g < 0) 
the homogeneous GPE supports bright soliton solutions. 



Vb(x, t) =y^sech(^^(a: ~ t;i))e^(— *)/^ 
^ (13) 
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where uj = ^mv + /z, /i being the chemical potential 
[1]. The bright soliton propagates with velocity v, and 
contrary to the dark soliton, its shape does not depend 
on its velocity, but only on the number of atoms in the 
condensate. 

In the dimensionless coordinates, the soliton solution. 



Eq. ([13|), is 

ipB{x,t) ^ Asech[B{x - vt)]e''^'"'-'^^\ (14) 

where A — B = N/2, the chemical potential of the con- 
^ , and w = ^ 



densate is /i = — and oj 
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The energy of the bright soliton, in the absence of an 
external potential is given by three contributions, '10]: 

^kin+-^int- (15) 
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The first part is the quantum pressure contribution to 
the kinetic energy 
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This is present even when the soliton is stationary, 
while the contribution to kinetic energy from the phase 
gradient 
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disappears for a stationary soliton. 

Finally, the interactions between the atoms contribute 

1 f IV"^ 
Eint^-^ dx I VI' 
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to the energy of the soliton. 

In total the energy of the bright soliton is 
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C. Dark-Bright solitons 

A dark-bright soliton is a composite object in a two- 
component condensates consisting of a bright soliton 
BEC and a density depression in a repulsive BEC. Here, 
we consider the case where all interactions, internal and 
between the two components, are repulsive with equal 
strengths, and in the rescaled coordinates, the coupled 
Gross-Pitaevskii Equations of the system read 
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This non-linear, coupled set of differential equations has 
a dark-bright soliton solution: 



^p-Dix,t) — l^iY^sino: + Y^cosatanh [^(a; — a;o)]^e 
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(21) 
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with V = K tan a [TTl [T2] . 

The bright soliton can only exists in the repulsive con- 
densate because it is trapped inside the dark soliton. The 
dark-bright soliton is wider than the single dark soliton, 
and the larger the number of particles in the bright com- 
ponent, the wider it is due to the repulsive interaction 
between the two components. 

The energy of the two-component condensate can be 
separated into three parts. 
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where the energy of the dark soliton E-^, is defined by 
subtracting the energy of the dark background conden- 
sate. 
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i?B is the energy of the bright soliton 
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and -Eint is the energy caused by the interaction between 
the two components 



iJint = / da: IV'bI'IVdI' (27) 



= iVBAisin a -I- -A^BAiCos a. 
The total vector soliton energy can then be written as 

m 

Ebb = ^At^ + ^Nbk^ + ^Nbv^ + MiVg. (28) 
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IV. SCATTERING OF BRIGHT SOLITONS 

Bright solitons are localized objects and in weakly 
varying potentials where the non-linearity is important, 
they largely behave as classical particles [T3HI8] , whereas 
in strongly varying potentials we expect wave diffraction 
phenomena to enter [101 ITSl - BS] . In this Section, we shall 
investigate the transition between particle and wave be- 
havior of bright solitons scattering on simple square po- 
tentials in one dimension. 



A. Weak Barriers 

As mentioned in the previous section, the chemical po- 
tential of a bright soliton is given by /i = — =^ , and when 
the potential is weak compared to the chemical potential, 
N"^ ^ Vex, we expect the strong interaction potential to 
hold the atoms together throughout the scattering pro- 
cess. 

In figure [l] the time evolution of the Gross-Pitaevskii 
wave function is shown for soliton wave packets, incident 
on a weak square barrier with four different incident ve- 
locities. The calculations show that the soliton keeps its 
sech-shape and is either reflected or transmitted by the 
barrier like a classical particle. 

In Figs. |lb| and[Tc|the soliton is slowed down dramat- 
ically and dwells on the potential for a rather long time, 
before it is eventually reflected or transmitted regaining 
its initial shape and speed. 

We define reflection and transmission coefficients as the 
weight of the wave function on the left and right hand 
side of the barrier 
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where L is the length of the box used for our calculations, 
and the barrier is centred dX x — 0. In Fig. [2] transmis- 
sion and reflection coefficients are shown for two different 
barrier widths. The figure confirms that the solitons are 
almost exclusively reflected or transmitted by the wide 
barrier, while the narrow barrier causes a minor splitting 
of the mean fleld wave packet. 

To get a better understanding of the observed phe- 
nomena, we assume that the soliton keeps its sech-shape 
during the scattering process, a fair assumption as seen 
from the numerical calculations in Fig. [T] 



ip{x, t) = Asech [B{x - xo)] e**«(^-^")+'''^, 



(30) 



and we use the Lagrange approach presented in |14l - 
116] to identify equations of motion for the time dependent 
parameters A, _B, xq and (j). The effective Lagrangian of 



the system is 
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= ^Nil -N<P+ ^NB{N -B)- V,s{xq), 
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(31) 



(32) 



We are now able to find the classical equations of motion 
for the soliton from the Euler-Lagrange equations 



dt dq dq ' 



q = B,xo, 



(33) 



Eq.(33) with q = yields iV = 0, i.e., the particle number 



is conserved. 



Eq. ( 33 ) with q — xq yields the equation of motion for 



the location of the soliton 



Xq = 



1 dV^sixo) 
' N dxo 
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This is the equation of motion of a classical particle with 
mass N and center of mass Xq moving in the effective 
potential, VcS, defined by Eq. (32 1. 



Finally Eq. (|33| with q = B yields 

,dVcs{xo) 
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describing a variation of the soliton width as it propa- 
gates within the potential barrier. 

For a square barrier, i.e., an external potential of the 
form 



Vix) 



for |a;| > § 
Vo for|x|<|. 
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the effective potential becomes 



VMxo)=< ^(tanh [B{l-xo)] 



for \x\ > I 



(37) 



-tanh -f xo)]) for |a;| < 



The Hamiltonian of the motion is found from the La- 
grangian: 

H = Nil-C= ^Nxl + ^NB{B - N) + V,s{xo), (38) 

and the threshold incident energy the soliton needs to 
cross the barrier is given by 



1 2 ^ 

-Nv? 

2 * 24 



max 



\xo\<{ 



NB{B -N) + Kff (xo) 



(39) 



5 



-50 



-25 





X 



400 



25 



50 



-50 



-25 



"T" 


X 



^400 



:^oo 



25 



50 




(c) 



(d) 



Figure 1: (color online) Time evolution of the condensate density. The external potential, indicated by the dashed 
lines, is a square barrier with height Vo = 0.04 and width 6 = 10 and the incoming soliton has = 1. The threshold 
velocity, Eq. (39), is vt « 0.289. In (a) the incoming soliton velocity is v 
,he velocity is v 
are v 



velocity, in|(b)|t^ 



0.24, which is well below the threshold 
0.29, which is just about the threshold velocity and in (c) and |(d)| the velocities 
0.3 and v 



0.4 respectively, which are above the threshold velocity. 



The corresponding threshold velocity of the bright soliton 
is indicated by the dashed lines in Fig. [2j and reproduces 
the transition between complete reflection and transmis- 
sion by the potential. 



We also note that since the barrier is softened by the 
interparticle interactions, Eq. (37), it makes sense that 



the soliton penetrates into the classical forbidden areas 
as seen in Fig. [l] [13]. 



B. Strong Barriers 



Scattering of bright solitons on strong barriers with 
Vex ^ A* tia-s been studied in the limit of high soliton 
velocities [lOl [HI [22l [23]. In this regime the soliton 
wave function may coherently split into reflected and 
transmitted components, and the asymptotic behavior 
of the transmission and reflection as the soliton velocity 
increases is the same as for a single particle, described by 
the linear Schrodinger equation. Here, we shall consider 
scattering of bright solitons on strong potentials with spe- 
cial attention to the behavior at low soliton velocities. 

Fig. [3] shows the time evolution of a bright soliton 
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Figure 2: (Color online) Transmission and reflection coefficients of a soliton scattering on a weak square barrier for 

various soliton velocities found by a numerical solution of the GPE (symbols) and by the Lagrangian method 
(dashed lines). The soliton has iV = 1 and the barrier h eigh t is Vq = 0.04. In (a) the barrier has width b = 10, and is 
thus wide compared to the soliton width = 2. In (b) the barrier has width 6 = 2 which is comparable to the 



soliton width. 



incident on a rectangular potential barrier. We observe 
both a reflected and a transmitted component in the cal- 
culation. When the potential barrier is high enough to 
split the soliton, i.e. Vq S> /i, and at the same time 
wide compared to the soliton width, a part of the soliton 
is seen to be transiently trapped inside the barrier, see 
Fig. [4j This trapping inside a potential barrier is due 
to the wave packet nature of the soliton. The trapped 
component oscillates and is gradually emitted out of the 
barrier region. This behavior is not seen when the barrier 
is narrow compared to the soliton width. 

We have solved the GPE numerically for different 
cases, and Fig. [5] summarizes our results (symbols) as 
functions of the incident velocity. When the external po- 
tential is strong compared to the chemical potential of 
the soliton, we expect the interaction term in the GPE 
to be less important, and the results in Fig. [5] may be 
understood from a linear approximation based on the 
Schrodinger equation. The solution to the Schrodinger 
equation for a single particle with energy E impacting a 
square potential barrier as in Eq. ( 36 ) is known, and the 
transmission coefficient is given by 



T =e{E)x 



iE{Vo~E) 



4E(E-Vo) 



sinh 



for E <Vo (40) 
sinh^ (^by^2{E - Va)^ ^ 
for E > Vo, 



where the Heaviside step function, Q{E), expresses that 
the energy of the particle scattered must be positive, 
which is evident in the single particle case, but will be of 
importance in the following. Ignoring the non-linear in- 
teraction term in the GPE is not valid in regions of space 
where the strong potential is not present. In these regions 
the soliton behavior is well known, and in particular we 
know the average energy per particle e: 
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(41) 



Because of the attractive interactions between the parti- 
cles, e becomes negative when v < (we have chosen 

= 1 in our calculations which gives « 0.29 in our 
dimcnsionless units). 

To describe scattering of the soliton in the linear ap- 
proximation, it makes good sense to replace the energy 



E in Eq. ( 40 ) b y th e average energy per particle in the 
' ' ' jV 

^/l2 



soliton, e in Eq. (41 ), and in the velocity regime v < 
the Heaviside stepfunction in Eq. 



(40 1 suppresses the 



transmission and predicts full reflection of the soliton. 
The prediction based on this description is shown in Fig. 
[5] and is in very good agreement with the solution of the 
GPE. The small transmission probability observed be- 
low the velocity threshold may be due, in parts, to our 
use of wave packets with finite spread, and hence with a 
momentum distribution exceeding the threshold value. 

We expect the energy per particle of the reflected and 
of the transmitted part of the soliton to equal the energy 
per particle of the initial soliton, ej. Since the transmitted 
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Figure 3: (Color online) Time evolution of the 
condensate density. The external potential is a square 

barrier placed at a: = 0. The potential has height 
Vo = 0.8 and width b = 0.5, and the incoming soliton 
has V = 0.5 and = 1. 
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Figure 4: (Color online) Time evolution of the 
condensate wave function. The external potential is a 

square barrier placed at a: = 0. The potential has 
Vq = 0.2 and b — 10, and the incoming soliton velocity 
is V = 0.7. Here we have zoomed in on the time 
evolution of the wave function |'0(2;)| (not the 
density), in order to see the details better. 



(reflected) parts contains a fraction Nx = T ■ N {Nj^ = 
R ■ N) of the atoms, the interaction energy per particle 
acquires a new value, and we thus expect the velocity of 
the components to change according to the equation 



coupled GP equations for this problem are then given by 
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Fig. [6] shows examples of the energy per particle in the 
transmitted and reflected components, determined by the 
GPE and by the expression ( 42 1 . There is good agree- 



ment between the energy per particle in the incoming and 
in the strongest of the outgoing components, i.e., with 
the reflected (transmitted) component for low (high) ve- 
locities. The poor agreement for the weakest outgoing 
component is a signature that it is not a pure soliton 
solution. 



V. SCATTERING OF DARK-BRIGHT 
SOLITONS 



Dark and dark-bright solitons retain their localized na- 
ture when they collide or propagate in slowly varying 
potentials [Tl[nill2, 26 30 . 

We will here study collisions of dark-bright vector soli- 
tons on the rapidly varying square potential, which we 
choose to only directly affect the bright component. The 
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+ \^uix,t)\^ + V{x) ]iPBix,t)i^B{x,t). 



Because of the rapid spatial variation at the edge of the 
square potential, we expect the vector soliton to break 
up into a reflected and a transmitted part as we observed 
for the single bright soliton. Due to the interaction with 
the dark component, this system is more complicated, 
and the linear approximation successfully applied to the 
bright soliton is not valid here. 

The results, shown in Fig. [7j are obtained by numerical 
simulations of the system (in all numerical calculations 
fj, = 1). The transmission and reflection coefficients, Tb 
and i?B, are defined in the same way as for the single 
bright soliton in Eq. (29 1. For the dark component. 



we define the reflection and transmission coefficients as 
the number of atoms missing with respect to the homo- 
geneous background condensate, on the reflection and 
transmission side of the barrier, divided by the missing 
number of atoms Njj in the incident dark component of 
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Figure 5: (Color online) Transmission and reflection coefficients of a soliton with TV = 1 on a strong barrier for 
vari ous incoming soliton velocities. The coefficients has bee n fo und by solving the GPE numerically (symbols), by 

(dashed lines), see the text. In 



Eq. (40) with E 



1,,2 



(dashed-dotted lines), and from Eq. (401 with E 



(a) flie barrier has strength Vq = 0.6 and width b — 2, in|(b)|the barrier has strength Vq — 0.2 and strength b — 2 



and in 



24 



(c) the barrier has Vb — 0. 8 an d b = 0.5, and thus the same area as in Fig. |(b)[ All barriers are higher than 



the chemical potential, = i. In (d) the barrier is a delta barrier of strength a = 0.4 which means that the barrier 



has the same area as the barriers in |5b| and [Sq For the delta potential the transmission probability reads 



T = e{E) 



/2E 



the vector soliton. 




The limits Xsw.t and Xgw.R are determined such that the 



amplitude variations due to sound waves emitted into the 
■^D-component during the scattering process do not con- 
tribute. We note from Fig. |8]that both the bright and 
the dark soliton are indeed split by the potential even 
though only the bright component is directly affected by 
the barrier. The dark soliton, however, is not reflected 
to the same degree as the bright soliton, and when the 
bright soliton reflection reaches unity, the dark soliton 
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Figure 6: (Color online) The energy per particle as a 
function of the incoming soliton velocity. After the 
scattering the energy per particle is found by 
numerical simulations of the GPE and Eq. ( 42 ) 



(symbols) and the energy per particle of the incoming 
soliton (solid line) is found by Eq. (41 1. 



transmission converges to a constant which is indepen- 
dent of a further increase in the barrier height, see Fig. 
[Sj This is also the case for other properties of the sys- 
tem such as the velocities of the outgoing solitons and 
the amplitudes of the sound waves emitted in the dark 
component. 

In Fig. [8] it is also seen that the reflection coefficient 
and the transmission coefficient of the dark soliton do 
not add up to one when the barrier height is increased. 
This we ascribe to sound waves carrying a corresponding 
surplus of atoms in the background condensate. 

When Nb is changed, it alters the soliton scattering 
properties since the coupling to the bright soliton, me- 
diating the effects of the barrier, is stronger when A'b is 
larger. 



rapidly at the edges, and where a more quantum mechan- 
ical behavior of the soliton is observed. 

For the bright soliton we found that the scattering 
properties of the soliton on a localized potential resem- 
bles that of classical or quantum mechanical scattering 
depending on the height and width of the barrier. 

When the barrier is low compared to the chemical po- 
tential of the soliton, the soliton behaves like a robust 
classical particle, since the barrier cannot overcome the 
internal interactions of the soliton. The potential ap- 
pearing in the resulting classical equation of motion de- 
scribing the soliton center of mass motion is an effective 
potential modified by the non-linearity. 

When the barrier is high compared to the chemical 
potential of the bright soliton, the soliton behaves like a 
single quantum mechanical particle with a modified rela- 
tion between its energy and velocity. In this regime, the 
internal interactions are no longer able to keep the soliton 
intact and it is split into a refiected and a transmitted 
part. 

The scattering properties of a dark-bright soliton on 
a barrier only affecting the bright component could be 
expected to resemble that of a single bright soliton scat- 
tering on a barrier. But we find that the scattering of 
a dark-bright soliton is in fact very different. The dif- 
ferences is due to the interaction between the two com- 
ponents of the dark-bright condensate, which enables an 
indirect interaction between the dark soliton and the ex- 
ternal potential. This interaction is strongest when the 
bright component is large compared to the dark one. 

Experimental realization of solitons in ID-condensates 
is well developed ^] , and dark-bright solitons have pre- 
viously been studied in slowly varying potentials. Thus 
the prospects of experimental studies along the lines of 
our calculations seem quite feasible. 



VI. CONCLUSION 



We have investigated scattering of solitons in Bose- 
Einstein condensates on localized potentials. The homo- 
geneous Gross-Pitaevskii equation supports matter wave 
solitons, and we have focused on two types of solitons: 
Bright solitons, which are wave packets without disper- 
sion due to the attractive, non-linear interactions be- 
tween the atoms in the soliton, and dark-bright solitons 
in two separate condensate components. In the latter the 
bright soliton is stabilized by the presence of the dark 
soliton, a dip in the condensate density. 

In slowly varying potentials, solitons are known to stay 
intact and their center of mass motion can be described 
by a classical equation of motion. We have here consid- 
ered scattering of solitons on square barriers which vary 
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Figure 7: (Color online) Time evolution of the condensate density when the dark-bright soliton scatters on a barrier. 



(a) shows the bright component and |(b)| the dark component of the condensate. The incoming soliton has = 3 
and sin a = 0.5, and the barrier, which only directly affects the bright component, has height Vq — 0.2 and width 

6 = 2, its location is indicated by the dashed lines. 
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Figure 8: Transmission and reflection coefficients of the dark soliton and the bright soliton as a function of the 
barrier height Vq. In (a) A^b = 3 and in (b) A^b = 1- In both cases sin a = 0.5 and the barrier width is that of a 



dark- bright soliton with sin a = (k, '-, Eq. (22)), for A^b = 3 that is 6 = 2 and for Nb = 1 it is 6 ~ 1.3. 
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